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Abstract. - 

We discuss the effects of collisions on the expansion of a degenerate normal Fermi gas, 
following the sudden removal of the confining trap. Using a Boltzmann equation approach, we 
calculate the time dependence of the aspect ratio and the entropy increase of the expanding 
atomic cloud taking into account the collisional effects due to the deformation of the distribution 
function in momentum space. We find that in dilute gases the aspect ratio does not deviate 
significantly from the predictions of ballistic expansion. Conversely, if the trap is sufficiently 
elongated the thermal broadening of the density distribution due to the entropy increase can 
be sizeable, revealing that even at zero temperature collisions are effective in a Fermi gas. 



The expansion of an atomic gas following the sudden removal of the confining trap is known 
to provide valuable information on the state of the system and on the role of interactions [1]. 
In particular in the superfluid phase, where the macroscopic dynamics of the gas are governed 
by the equations of hydrodynamics, one expects anisotropic expansion if the gas is released 
from an anisotropic trap. This peculiar feature, which should be contrasted with the isotropic 
ballistic expansion exhibited by a non interacting gas, was first discussed in [2,3] in the case 
of Bose gases and in [4] in the case of Fermi gases. Anisotropy is not however a unique feature 
exhibited by superfiuids and also in the normal phase one can expect a similar behaviour if 
collisions are sufficiently important [3,5-7]. This effect is expected to be particularly important 
in the case of Fermi gases interacting with large scattering lengths near a Feshbach resonance 
[8-10]. At first sight one would expect that the effects of collisions are suppressed at low 
temperature because of Pauli blocking. This is certainly true if one works close to equilibrium, 
as happens in the study of small amplitude oscillations [11]. In the problem of the expansion, 
however, large deformations in momentum space can be produced if one starts from a highly 
deformed cloud, with the result that collisions become effective even if the gas is initially at 
zero temperature. This interesting possibility was first pointed out in [12]. 

The aim of the present work is to calculate explicitly the dynamics of the expansion of 
a dilute, degenerate normal Fermi gas, taking into account the role of collisions. The main 
purpose is to provide quantitative predictions for the aspect ratio and the thermal broadening 
of the density distribution, as a function of the relevant parameters like the ratio of the trap 
frequencies, the scattering length and the number of atoms. 
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The system we consider is a dilute two component Fermi gas. At low temperatures the 
collisions between two atoms of the same species are suppressed due to Pauli blocking, and 
only atoms of different species can collide. We assume that the two species have the same 
mass and density. The starting point is the Boltzmann equation 
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where for each component /(r, vi,t) is the distribution function in phase-space and f7 ex t( r ) = 
(m/2)(uj 2 _{x 2 + y 2 ) + u 2 z 2 ) is the external trap potential. We note that the potential is 
cylindrically symmetric, which is the form favoured by experiments. Moreover, in we 
have neglected the mean field interaction term [13] which, however, has a minor effect on the 
expansion of a dilute Fermi gas [4] . The collisional integral in the case of a dilute Fermi system 
reads 
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where in the low-energy limit, a — Ana 2 , with a the s-wave scattering length. The scattering 
length is assumed to be smaller than the average distance between atoms. 

The dynamics can be studied analytically using the method of the averages [14,15] which in- 
volves calculating moments of the type (x)(t) — J x( r : v )/( r i v i t) d 3 r d 3 v/ J f(r,v,t)d 3 rd 3 v. 
Applying this to JQ) gives a set of equations of the general form 
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where x is a general function which depends on r and v. Here we shall take moments of the 
form (rf ) , (r^) and {(vi — Ui) 2 ), where Ui(r) represents the velocity fields within the gas. If 
X is a conserved quantity during a collision (i.e. A% = \i + Xi ~ Xv ~ Xi' is zero) then the 
collisional term disappears. This is true for the moments (rf) and (r^) as well as for the 
quantity J2i(( v i ~ u i) 2 )- This method has been used to study collective excitations in [11]. 

To explicitly include the collisional term in the calculations, it is convenient to introduce 
the following parametrization for the distribution function 



/(r,v,i) 



exp 



£?(r,v,i) - n 
T 



(4) 



with 



E 



(Vj - Uj) 2 



(5) 



where Ui — PiTi, and Pi, bi, Ki, /i and T are time dependent parameters, with i € {x,y,z}. 
The ansatz ll 1151 contains more free parameters than necessary, so that we can, without any 
loss of generality, set Yli^i — (Tli^'O -1 ^ 2 - This fixes a unique definition of the effective 
temperature T [16]. The chemical potential, //, is then obtained from the normalization 
condition, (T /fiuj) 3 f^[e^^ T ) = N/2, where N is the total number of atoms in both components, 
d> 3 = Jli^j' an d fs( z ) = J2 n (~ l)" +1 z"/?7 ,s . From Eas. (|4l5|) one can evaluate the entropy per 
particle of the gas, S = (4p-/x)/T, where p = T/ 4 (e M/T )// 3 (e Al/T ). This result, taken together 
with the normalization condition, shows that there is a one-to-one correspondence between 
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the entropy and the temperature T, and consequently changes in T taking place during the 
expansion are associated with a change of entropy. 

The ansatz (1 1151) includes the initial equilibrium configuration predicted by Fermi statistics, 
and accounts for rescaling effects in coordinate as well as in momentum space. It allows, in 
particular, for anisotropic effects in momentum space which are crucial for describing correctly 
the mechanism of the expansion. The parametrization can describe both the collisionless and 
the hydrodynamic expansion as well as other intermediate regimes. Furthermore it accounts 
for possible changes in the effective temperature of the system during the expansion. 

The next step is to evaluate the moments (x) in terms of the scaling parameters, which 
gives (r?) = bfp/muj 2 . Other moments can be expressed in terms of this quantity, so that 
(riVi) = (3i{rf), ((v l - « t ) 2 ) = Lo 2 Ki{r 2 )/b 2 , and {u 2 } = f3f(r 2 }. The moment equations then 
become 

A+/f-|Mc+«?-ft (7) 
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where the final term in Q) arises from the confining potential. We set this term to zero in 
studying the expansion of the gas, while rescaling time in units of u>±, the original radial 
trap frequency. In addition, it is convenient to rewrite © as equations for the anisotropy in 
momentum space, s = (K z j 'ifjj 1 / 2 , and the entropy, S 
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where we have defined 
with J(s,r) given by 
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Here V — (V\ x — V2 X ) 2 + {V\ v — Viy) 2 + s 2 (V\ z — V~2 Z ) 2 , and we have rescaled variables so 
that /(V) = [exp(i? 2 + V 2 — /i/T) + l] -1 . The function J depends on both s and the reduced 
temperature t — T/Tf, where Tf = (SN) 1 ' 3 ?^!) is the Fermi temperature. In Eqs. 1|9I11J) we 
have introduced the relevant dimcnsionless interaction parameter 

£ = (AA0 1/3 (M 2 , (13) 

where kp is the initial Fermi wavevector at the center of the trap, and A = lu z /luj_ is the trap 
anisotropy. Further, we have introduced the geometric average K = (K 2 K z ) 1 l' i = i\) 2 x b z )~ 2 l z . 

Eq. IjlUI) explicitly shows that the entropy (and hence the temperature) of the gas will 
remain constant during the expansion (S = 0) either in the absence of collisions (J = 0) 
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or when the distribution in momentum space is isotropic (s = 1) [17]. The latter situation 
arises when starting from a spherical trap, A = 1, or when collisions are sufficiently frequent 
so that the gas is in the hydrodynamic regime. In regimes intermediate between the colli- 
sionless and hydroynamic limits, deformations in momentum space produce an increase of 
entropy and temperature. The equations (l(ill2ll can also be used to study small oscillations 
around equilibrium, where the momentum distribution is isotropic (s = 1). This proce- 
dure reproduces exactly the results of [11], in particular the collisional term takes the form 
((^- Ul ) 2 C//}/((^- Ml ) 2 } = («-!)/?, with 1/(wit) = CK 2 £F q (T/T f ) andC = 8/(5 -3 5 / 3 ). 
Here f is the typical relaxation time for the shape oscillations and Fq is given in [11]. Note 
that 1/f a T 2 at small temperatures, so that the collisional contribution disappears when 
T = 0. This is the result of Pauli blocking of collisions for a spherical momentum distribution. 
However, during expansion large deformations in momentum space can lead to collisions that 
scatter atoms outside of the Fermi surface, and the function J differs from zero, even at zero 
temperature, as pointed out in [12]. 

We will assume that the gas is initially in equilibrium at low temperature, and that the 
effective temperature contained in the ansatz Q remains small during the expansion. One 
can then expand the Fermi functions f s (z) using low-T Sommerfeld expansions, so that the 
normalization condition gives \x ~ T F (1 — 7t 2 t 2 /3), p ~ T F (l + 27r 2 r 2 /3)/4 and S = tt 2 t. 
Substituting these expressions into (|6I7I9I10[) yields a new set of equations that can be solved 
numerically to study the time evolution of the gas. To simplify matters we also approximate 
the function J(s, r) with its zero temperature value 

2 • 3 4 / 3 s 5 / 3 r 

J M) = ~6 / d 3 Rd 3 V 1 d 3 V 2 d 2 nV'AV z 2 9o(V 1 )e (V 2 )(l - 0b(Vi))(l - 0o(V£)), 

where we have rescaled the coordinates again such that 9q(V) = Q(l — V 2 — R 2 ), where O is the 
Heaviside step function. The integral for J(s) can be evaluated numerically using a standard 
Monte Carlo technique, and is plotted in Fig.^ Also plotted is the result J(s) = 2s 8//3 /(3 5 / 3 7r) 
holding for large s, where the integral can be evaluated analytically by noting that virtually 
all collisions in this limit will scatter atoms outside of the Fermi surface (i.e. (1 — ^o(V' 1 ))(l — 
flo (V2)) — 1). We shall discuss the validity of this zero-temperature approximation later. 

We now use these equations to study the expansion for different choices of the trap 
anisotropics, A < 1, starting from a T = configuration [18]. We plot the aspect ratio 
of the cloud \J (rj) / '(z 2 ) = Xb± /b z as a function of time for £ = 1 and £ = 20 and compare to 
the expected collisionless and hydrodynamic behavior in Fig. [21 One notices that the effects of 
collisions are far less pronounced for A = 0.3, where the results are barely distinguishable from 
the collisionless behavior. However, for A = 0.03 the expansion lies intermediate between the 
two limits if one chooses £ = 20. This is to be expected since the deformation in momentum 
space reached in this case is much larger. In particular, in the collisionless case one would 
expect s — ► 1/A at long times, and this fixes the maximum s possible. Since J takes large 
values only for large s, then the collisional term is potentially more effective for smaller A. 

An interesting consequence of the momentum space deformation is that atoms scattering 
outside the Fermi surface will tend to smooth out the distribution function. This aspect 
appears in (|10fl as an increase in the temperature, even if we start from an initial T = 
configuration. Fig. Ofa) shows the temperature calculated at uo\t — 100 for different values 
of A. Again, one sees a much larger effect for A <C 1. The function tends to zero for both 
£ — > and £ — > 00, representing the crossover from collisionless to hydrodynamic behavior in 
a similar manner to damping times in collective oscillations [11,14]. Since l|10|) is an equation 
for t 2 one finds that starting from a non-zero temperature, r , gives a final temperature of 
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approximately r = -\/tq + {St) 2 , where St is the change for an initially zero temperature gas 
obtained from Fig. OJa) . Hence the temperature increase produced by collisions will become 
smaller as one raises the initial temperature of the gas. 

The change in temperature during the expansion, as well as the fact that experiments will 
be initially at a non-zero temperature, will lead to corrections to the zero-temperature colli- 
sional term (| 1 4|> . We can estimate the effects at hnite temperatures by approximating J(s, r) 
Ijl2|l with the sum of the zero-T result to that derived for finite-T but small deformations, so 
that J(s, r) = J(s, 0) +C(s - I)Fq(t), with the function Fq(t) given in [11]. We find that the 
inclusion of the temperature dependence in J has little impact on the results, even allowing 
for an initial temperature of r = 0.2. 

Since for a dilute gases (kp \a\ <C 1) the value of £ for realistic parameters will be, at most, of 
the order of 1 , we conclude that collisional effects on the aspect ratio of the expanding gas are 
negligible (see Fig. I2J • This is consistent with experimental results in dilute degenerate gases 
(see for example [9,19]). In contrast, we find that the effect of collisions on the thermal broad- 
ening can be significant at low temperatures even if £ ~ 1, especially for highly deformed traps. 
As an example in Fig. Efb) we show the column density n(r±) = (m/h) 3 J dzdvf(r,v,t), 
calculated at uj±t = 100, A = 0.03 and £ = 1, starting from an initial zero temperature 
configuration. Around 80% of the increase in T (and the consequent broadening of the den- 
sity profile) takes place over the first LJ±t = 20 of the expansion. The comparison with the 
prediction of ballistic expansion (dashed line) explicitly reveals the importance of collisions. 
Experimentally one could observe this difference by cooling down a two component Fermi gas 
to very low temperatures. Ballistic expansion could be achieved either by first removing one 
of the two components, or by suddenly tuning the scattering length to zero at the start of the 
expansion. 

In order to also observe large effects in the aspect ratio one should increase the value of 
£, and hence of fecial, by, for example working close to a Feshbach resonance. In this case, 
however, the formalism of the Bolzmann equation is not strictly applicable. A rough estimate 
of the collisional effects can be obtained by replacing a 2 with the unitarity limited expression 
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Fig. 2 - Aspect ratio against time for the free expansion of a zero-temperature Fermi gas from a 
trap with (a) A = 0.3 and (b) A = 0.03. The upper and lower solid lines on each plot represents the 
behavior in the hydrodynamic and collisionless regimes respectively, with the other lines displaying 
results of solving the equations 161121 for f = 1 (dashed) and £ = 20 (dotted). 



a 2 / (1 + K (kpa) 2 ), where Icf is the initial Fermi momentum and K accounts for the decrease 
of the density during the expansion. In the unitarity limit kp\a\ — > oo equation is then 
modified by replacing K 2 with K and setting £ = (XN) 1 ^ 3 . These changes result in a sizable 
increase of the anisotropy effects. Apart from the fact that the parameter £ can easily take 
large values, the replacement of K 2 with K makes the collisional term effective for longer times 
during the expansion. One should however note that in the unitarity limit the gas is expected 
to be superfluid at low temperatures and its dynamics should be consequently described by 
the hydrodynamic equations of superfluids. 

In conclusion we have shown that collisions can be effective in a dilute normal Fermi gas 
even at zero temperature, as a consequence of large deformations of the distribution function 
in momentum space after expansion from a very elongated trap. They can give rise to a 
sizeable entropy increase and hence thermal broadening of the density distribution, which 
should be visible by imaging the atomic cloud. In contrast a T = superfluid should expand 
anisotropically, without any entropy increase due to the absence of collisions. 
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